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A dynamical system framework is used to describe transport processes in plasmas 
embedded in a magnetic field. For periodic systems with one degree of freedom the 
Poincare map provides a splitting of the phase space into regions where particles have 
different kinds of motion: periodic, quasi-periodic or chaotic. The boundaries of these 
regions are transport barriers; i.e., a trajectory cannot cross such boundaries during 
the whole evolution of the system. Lagrangian Coherent Structure (LCS) generalize this 
method to systems with the most general time dependence, splitting the phase space into 
regions with different qualitative behaviours. This leads to the definition of finite-time 
transport barriers, i.e. trajectories cannot cross the barrier for a finite amount of time. 
This methodology can be used to identify fast recirculating regions in the dynamical 
system and to characterize the transport between them. 


1. Introduction 

Transport phenomena are ubiquitous in nature and involve the redistribution of 
physical quantities such as mass, charge, linear and angular momentum and energy etc. 
Different mechanisms are at work in the transport processes ranging from diffusion to 
advection and mixing in the case of turbulent or chaotic motions. Recently a new concept 
in the study of transport processes in complex fluid flows was introduced by Peacock & 
Haller (2013): Lagrangian coherent structures (LCS). In a two-dimensional configuration 
these structures are special linesf advected by the fluid which organize the flow transport 
processes by attracting or repelling the nearby fluid elements over a finite time span. 
Strictly speaking we should call these structures Hyperbolic LCS. Two other kinds of 
LCS have been introduced in the literature. For the definition of all the different kinds 
of LCS we indicate the recent review by Haller (2015). For the sake of simplicity we will 
refer to these structures simply as LCS. These special lines organise the flow splitting the 
domain into macro-regions with fast mixing phenomena inside them. Over the finite time 
span which characterizes the LCS these macro regions do not exchange fluid elements. 
The LCS have been widely used in the literature to characterize transport processes in 
various systems: the pollutant transport on the ocean surface (Coulliette et al. 2007), 
blood flow (Shadden & Taylor 2008), the spreading of plankton blooms (Huhn et al. 
2012), turbulent combustion (Hamlington et al. 2011), jellyfish predator-prey interaction 

f Email address for correspondence: falessi@fis.uniroma3 

f In the most general system the LCS are surfaces advected by the flow but here for the sake 
of simplicity we will only deal with 2D systems. 
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(Peng & Dabiri 2009), atmospheric dataset analysis (Tang et al. 2010), solar photospheric 
flows (Chian et al. 2014), saturation of a nonlinear dynamo (Rempel et al. 2013), etc. 

Plasmas are often studied using “fluid” theories either in phase space, such as the 
Vlasov-Maxwell system, or in physical space such as the two fluid and the MHD systems. 
The LCS techniques can therefore be applied to study transport processes, i.e. the mixing 
of fluid elements, in these systems. In (citation to another proceeding: Carlevaro, Falessi, 
Montani, Zonca) the LCS has been used to quantify the phase space transport due to 
the interactions between two supra-thermal electron beams and a cold, homogeneous, 
background plasma. In two recent works (Borgogno et al. 201 la, 6 ) it was shown how the 
LCS can provide information about the electron transport due to the stochastization of 
the magnetic field in a collisionless reconnection process. 

The introduction of these techniques is relatively recent and, in spite of their increasing 
use, their rigorous definition has been subject to debate (Shadden et al. 2005; Haller 
2011). The first definitions of an LCS were based on the Finite time Lyapunov exponent 
profile. These have been shown to be incorrect by G.Haller who found several coun¬ 
terexamples to this heuristic definition. The rigorous definition of an LCS as the most 
repulsive or attractive material line with respect to the nearby ones was introduced by 
(Haller 2011). In this article we provide a simplified version of this derivation in order to 
give to the reader some physical intuition about these structures, we analytically calculate 
the shape of the LCS in a simple Hamiltonian system and, finally, we compare the LCS 
obtained by Borgogno et al. (20IIa) using the FTLE method with the ones obtained 
with the rigorous definition introduced. The numerical framework used to compute the 
LCS shape was developed by Onu et al. (2015). 


2. Lagrangian Coherent Structures and Transport Barriers 

2.1. Lagrangian Coherent Structures as most repelling material lines 
Following (Haller 2011), in this section we briefly review some mathematical concepts 
that lead to the definition of LCS. 

We consider a dynamical system in 2D phase space x = {x,y), 

dx , . dy , . , . 

— = u„(t, x,y), — = Vy{t, X, y) (2.1) 

with continuous differentiable flow map 

= x(t,to,a^o)- (2.2) 

Two neighbouring points xq and xg -I- Sxg evolve into the points x and x -|- dx under the 
linearized map 

|fe) = V< |feo). (2.3) 

where for notational convenience, we adopt a bra-ket notation for vectors and scalar 
products and represent a generic column vector as |c) and a row vector as (r|. Their 
scalar product is denoted as d = (c|r). 

Consider a curve 70 = {xq = 'r(s)}. At each point a:o G 70 : define the unit tangent 
vector |eo) and the normal vector |no). In the time interval [tojf] the dynamics of the 
system advects the “material line” 70 into 74 and Xq G 70 into Xt £ ■ The linearized 

dynamics maps the tangent vector |eo) into \et) which is tangent to jt 


let) 


'^(l>loixo)\eo) 


eo 




eo 


V</»L(a;o)|eo) 
^(eo |C'to(a^o)| eo) 


(2.4) 
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where Cl^{xo) = ^'Ato Cauchy-Green strain tensor and ^ stands for 

transposed. This symmetric tensor describes the deformation of an arbitrarily small circle 
of initial conditions, centered in Xq caused by the flow in a time interval As an 

example let us consider a circle centered in a;o with radius ||5a:o||- After the time interval 
[to, t] this will be deformed into an ellipse with major axis in the direction of \^rnax) and 
minor axis in the direction of |^mm) being \^rnax) and |^mm) the two eigenvectors of 
Cfg(xo). The corresponding real and positive eigenvalues are Amax and Xmin- The length 
of the major axis is Amoz,||<5iCo|| while the the length of the minor is Aminll^a^oll- The 
curves with tangent vector along |^mm) and, respectively, \^rnax) are called strain lines 
of the Cauchy-Green tensor. In general, the mapping will not preserve the angle between 
vectors and therefore usually \nt) differs from ^4>lg |no). 

Using the condition (no|eo) = (no | | eo)=0 and inserting Eq. (2.4) we obtain 

the expression for \nt) which is given by 


Wt) 


{Vct/°Y\no) 


no 




no 


(V</>[°)^|no) 
A/(no |C'-i(a:o)|no)’ 


(2.5) 


where C~^{xo) = Cf°(xo) and the time interval marks have been suppressed as will be 
the case in the following formulae when not explicitly needed. We define the repulsion 
ratio p\^(xq, no) as the ratio at which material points, in other words points advected by 
the flow, initially taken near the point a^o G 7 o, increase their distance from the curve in 
the time interval [to,t]: 

pl(a^o,no) = (nt |V 0 *^(a;o)| no). ( 2 . 6 ) 

In order to understand this definition we may imagine measuring the distance between 
Xq, taken on the curve 70 , and a point initially placed at unit distance from the curve. 
After a sufficiently small amount of time this distance is measured by the projection of 
the vector V0)g(a:o)|no) along \nt) as shown in Figure I. Using the previous definitions, 
pl^{xo,no) can be expressed either in terms of no or of n* as 

pl{xo,no)= \ = A/(nt |C(a;o)|nt). (2.7) 

a/( no \C-^{xo)\ no) 

Similarly, the contraction rate L\^(xq) is proportional to the growth in time of the vector 
tangent to the material line 

L{xo,eo) = 1/(60 |C'(a:o)| eo). (2.8) 


The aim in these definitions is to characterize a LCS over a finite time interval [to, to + T] 
as a material line, in other words a curve advected by the flow, which is locally the 
strongest repelling or attracting curve with respect to the nearby ones. This leads, as 
shown by Haller (2011), to the following definitions. 

Definition I. A material line satisfying the following conditions at each point: 

(i) 

^min ^max^ ^max ^ 1 (^*^) 

(ii) the tangent vector is along the eigenvector associated with the smallest eigenvalue 


l^o) — \^min) 5 


( 2 . 10 ) 
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Figure 1. The green line and the purple line represent the trajectories of two elements of the 
system. The continuous thick black line is the curve 70 while the dashed one is its evolution 74 
at time t. We have marked in blue the evolution of |no) through the linearized dynamic and in 
red its projection on \nt). 


(iii) the gradient of the largest eigenvalue is along the curve 




VA. 


= 0 


is called a Weak Lagrangian coherent structure (WLCS). 


( 2 . 11 ) 


Definition 2. A WLCS which satisfies the following additional condition 
(i) at each point the relationship 


(e 


max 


fmax I 


< 0 


holds is called a Lagrangian coherent structure. 


( 2 . 12 ) 


A simplified illustration of the previous conditions can be given as follows. The first 
condition is obtained when requiring that at each point Xq of the material line the 
repulsion rate p{xQ,nQ) is larger than the contraction rate L{xQ,nQ) which represents 
the effect of the shear along the material line. At each point along a material line the 
tangent vector can be expressed in terms of the eigenvectors \fmin) and \fmax) of the 
Cauchy-Green strain tensor C: 

|eo) = \fmin) + /3 \fmax) > + /3^ = 1. (2-13) 


where a and jd represent the orientation of the material line and may be arbitrary 
functions of Xq. It follows that 


|lIo) — O; \frnax) (d |^mm) ■ 


(2.14) 
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We can therefore express the repulsion rate in terms of a and /3 


pl^{xo,no) 


1 

a/( no \C-^\no) 



(2.15) 


Now at each point Xq we maximize pl^ with respect to the direction of the tangent vector 
|eo) i.e., of the orientation of the chosen material line, with the constraint + /3^ = 1 . 
This yields a = 1, /3 = 0, so that 


pI 


0 ~ 


(2.16) 


which implies that the material line must be chosen to be the strain line oriented along 
the eigenvector |^min) as required by Eq. (2.10). At this stage we need to maximize the 
repulsion rate along the material line with respect to nearby material lines. Therefore we 
define 

^‘ 0 (^ 0 )=/ dsp{xois),nois)), (2.17) 

J Si 

which has the physical meaning of the repulsion rate integrated over a curve 79 . We 
consider a curve 7 ^ with points x^ G 7 ,; such that 


Xf^ = X(j + eh{s) |no). 

The first variation of Eq. (2.17) with respect to Eq. (2.18) gives 
SRloile)[h] = 5 [ ds p{x^{s),n^{s)) = 0, 

J Si 

which can be computed as 

Q /•S2 Q j'S 2 _ 

lim — / dsp{Xe{s),ne{s)) =\\vn— / ds^/Xmax (Xe(s),ne(s)) 

de de 


(2.18) 

(2.19) 

( 2 . 20 ) 


leading to 

= r ^f^iVX^axl^max) = 0 ( 2 . 21 ) 

J Si V '^max 

which vanishes if \^Xmax) is tangent to the material line. Requiring that this material 
line represents a maximum of the integrated repulsion rate we obtain 


i^^max ^ Xjjiax ^max^ ^ 0 


( 2 . 22 ) 


which is the condition that defines the locally most repelling LCSs. These structures are 
Lagrangian by definition and have no transport through them because they are material 
lines. 


2.2. Lagrangian Coherent Structures as second derivative ridges 
Second derivative ridges were defined by Shadden et al. (2005) in terms of the features 
of the Lyapunov exponent field a{xQ,tQ,t) that characterizes the rate of separation of 
close trajectories. 

The finite time Lyapunov exponent can be expressed in terms of the Cauchy-Green 
tensor eigenvalues as 


a{xo,to,t) 


1 

2\t-to\ 


fo Xjnaxi^O: ^0i ^) 


(2.23) 


Definition 3. Curves r{s) (not necessarily material lines) such that 
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(i) the tangent vector \r'{s))and Vcr along the curve are parallel, 

(ii) the normal unit vector \n) is such that along the curve for all {u\u) = 1 

(n iri n) = min {u 1171 m) < 0 (2.24) 

where S is the Ftessian matrix of the second derivatives of a with respect to Xq is called 
a second derivative ridge. 

A major difference between the two sets of definitions is that the most repelling LCS 
definition involves the eigenvectors and eigenvalues of the Cauchy-Green strain tensor, 
while the second condition in the definition of the second derivative ridge is governed by 
the eigenvectors and eigenvalues of the matrix S. 


3. Example: a Hamiltonian flow map 

In this section we will illustrate the definitions introduced above and the procedure 
needed in order to identify and characterize the Lagrangian coherent structures by direct 
construction in a simplified dynamical system that can be studied analytically. We find 
it convenient to consider a Hamiltonian system which ensures condition (2.9) through 
conservation of phase space volume. We consider the time independent Hamiltonian 



(3.1) 


from which we obtain 



(3.2) 


A contour plot of the Hamiltonian FI is given in Fig. 2. The trajectories obtained by 
integrating Eq. (3.2) coincide with the lines of constant FI. Integrating the system (3.2) 
we obtain 


x{t) = {4>l{xo,yo)) ^ ( 3 - 3 ) 

y{t) = {4>l{xo,yo))y=yo{l + ^"^ , (3.4) 


where xq = x(t = 0) and go = y{t = 0). From Eq. (3.3) we see that points with negative 
xq reach a; = — oo in a finite time while points with positive xg reach a: = 0 for t —>■ oo. 
In order to avoid this finite time singularity we restrict the domain to the semi-plane 
Xq > 0. 

The matrix V^q takes the form 


= 


1/B^ 0 

Btyo 


(3.5) 


where B = 1 + txgl2. The Cauchy Green strain tensor CQ{xo,yo) = (V^g)^ '^4>o is 




Btyo 

B^ 


(3.6) 


In order to identify its eigenvectors and corresponding eigenvalues it may be convenient 
to diagonalize it 


Riif) C^o i?^(^) = D. 


(3.7) 
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Figure 2. Contour plot of the Hamiltonian shown in Eq. (3.1). 


by means of the rotation R{ip) 


R{ip) = 


COS (f — sin ip 
sin(p cosif 


with (p given by 


I f 2Btyo 


and D{xo,yo) is the diagonal matrix 


D = 


'min 

0 


with \max — -^+5 -^mm — S-nd 

1 + + BH^y^ ± yJ[l + B» + BH^y^f -AB» 


A+ = 


2B4 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


The rotation angle ip{xo,yo) can be used to find the orientation of the strain lines. Note 
that <p(xo, yo = 0) = 0, as the strain tensor in Eq. (3.6) is diagonal on the positive semi 
axis t/o = 0 with |^min) tangent this axis i.e., in the xq direction. In fact from Eqs. (3.3) 
we see that if we take two points xqi and xo 2 on the positive yo = 0 semi-axis, with 
2^02 > 2 ;oi they stay on this axis and their distance d(t) decreases in time: 


d(t) 


_ m _ 

1 + t (X 02 + Xoi) /2 -I- t^X02X0l/4:' 


(3.12) 
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Figure 3. Plot of the eigenvectors of the Cauchy-Green strain tensor in the limit txo 2. 
\^'max) is marked in red while |^mm) in blue. 


as consistent with the fact that this yo = 0 axis is parallel to |Cmin) at all times. In order 
to show the strain lines of the Cauchy-Green tensor for a;o > 0 in Fig. 3 we plot the vector 
fields \imax) = 1 ) = |(sin(/ 3 , cos(/ 3 )) and \^min) = 0 ) = I cos (/?, - sin (/?)) 

in the limit of sufficiently long times such that teg 2 in which limit the rotation angle 
(/? reduces to: 

(f ^ (1/2) arctan[(xo?/o)/(yo ~ ^Jn/d)] = arctan (3.13) 

Xo 

and becomes independent of time. In the same limit the eigenvalues of the Cauchy-Green 
strain tensor become simply 


^max ^ t ^Xq -t- I^qXq^ /I 65 ^min — ^ 


and have a factorized time and space dependence. 

From (3.14) we can compute 'VXmax (which will be parallel to ’VXmin) 

^ ^0 f Xq -j- 2^0 
2 xo?/o 


^ Xmax — 


OC COS ip 


cos^ (fi 2 sin^ (fi 


2 cos (fi sin (fi 


(3.14) 


(3.15) 


A contour plot of Eq. (3.14) is shown in Figure 4. Recalling the sufficient and necessary 


condition in order to have a WLCS is 




V A„ 


= 0 we find the condition 


cosip sin(2sin^ (/? -|- 3 cos^(/j) = 0. For Xg ^ 0 (i.e (p ^ this condition implies 

sin(/? = 0. Thus the only repulsive WLCS is given by the positive yg = 0 semi-axis. 

It is easy to see that the condition (2.12) is not satisfied, i.e., this WLCS is not a 
maximum for the repulsive rate and thus is not an LCS. 

Referring now to the definition of the second derivative ridges, from 3.14 at j/g = 0 we 
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Figure 4. Contour plot of Xmax- 


find 

S{xo,yo = 0) = ^ ^ J 2 ^ (3.16) 

which is not negative definite and for which {n\E\n) is a maximum, not a minimum, 
contrary to the requirement (2.24) in the second derivative ridge definition. 

In this example the WLCS according to Haller’s definition that we found at j/o = 0 is 
not a second derivative ridge. 


4. Electron transport in a 3D reconnection process 

In this section we will deal with the barriers to the transport of electrons in a 3D 
magnetic reconnection processes in the presence of a fixed “guide field” component . A 
simplified doubly-periodic slab geometry is considered. Because of this double periodicity 
this model can be applied mainly to toroidal thermonuclear plasmas (Borgogno et al. 
2005; Avinash et al. 1998; Porcelli et al. 2002), and with some modification to space 
plasmas (Rappazzo & Parker 2013) and, more generally, to physical contexts where 
plasma is confined by a nearly uniform strong magnetic field with small perturbations. 
Under proper adiabaticity conditions, and in particular assuming that their Larmor 
radius is negligible, charged particles move along magnetic field lines. Therefore in this 
approximation the study of the motion of electrons can be referred to the study of the 
topology of the magnetic field lines. 

The magnetic field structure and evolution are described through the evolution of the 
magnetic flux function y, z, t) by the relationship 


B = Boe^ + X V^. 


(4.1) 
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As for all solenoidal fields in an odd dimensional space, at any fixed time t = t the 
magnetic field line equations can be written in the form of Hamilton’s equations. The 
function ip{x,y, z,t) plays the role of the Hamiltonian while the coordinate z plays the 
role of the time variable in this Hamiltonian system. Using a periodic geometry along z 
allows us to study the shape of magnetic field lines with dynamical systems techniques 
such as the Poincare map. The growth and interaction of different unstable modes will 
naturally lead to the formation of structures in the magnetic field such as magnetic 
islands and to chaotic behaviour of field lines. Assuming, as commonly done, that the 
time it takes electrons to complete a number of turns along z before the magnetic field 
is significantly changed, we can investigate the advection of particles in such a chaotic 
system at a fixed t = t using the LCS technique. This requires that the electron thermal 
velocity be sufficiently larger than the Alfven speed and allows us to highlight the finite 
time transport barriers of the system. Choosing the LCS characteristic time span |z — zo| 
properly, we obtain the boundaries of the regions where fast electron mixing is expected. 
After a time span of the order of the Alfven time we can expect that the shape of the 
magnetic field lines will have changed significantly and therefore we need to “refresh” the 
Hamiltonian y, z, t) and plot a new Poincare map for the magnetic field. The LCS of 
the Poincare map calculated with a sufficient number of iterations marks the boundaries 
of the regions where the electrons will mix. 

In this paper we will not deal with the numerical integration of the PDE governing 
the dynamics. We will instead study the shape of the magnetic field lines at fixed t 
extracting from the numerical simulation carried out by Borgogno et al. (2005). An 
analogue analysis has been carried out by Borgogno et al. (201 lo) where the LCS have 
been calculated using the definition based on the ridges of the Finite time Lyapunov 
exponent field. We choose the same values for the parameters defining the LCS in order 
to make a comparison with this work. 


4.1. The physical system 

The physical system studied is a Hamiltonian reconnection process in a dissipationless 
3D plasma immersed in a strong, uniform, externally imposed magnetic field. A slab 
geometry is used with a formal additional periodicity along x imposed for numerical 
convenience. The algorithm applied in the numerical simulation is detailed in Borgogno 
et al. (2008). The reconnection process develops in a static equilibrium configuration 
given by: 

i^eq{x) = Acos{x), (4.2) 

with A = 0.19. The integration domain is defined by —L^ < x < L^, —Ly < y < Ly, 
—Lz < z < Lz with Lx = TT, Ly = 27r and Lz = IOtt. The equilibrium is perturbed as 


F{x, y, z, t) = tpeq{x) + tp{x, y, z, t), (4.3) 

where if) is written as a sum over Fourier modes: 

ip{x, y, z, t) = Si 'tj)i{x, kyty + kzrZ, t) (4.4) 

with kyi = kzi = The field line equations can be cast in Hamiltonian form: 


dx _ 1 d'l' _ Bj; 

dz ~ Bq dy ~ Bq 

dy ^ 

dz Bq dx Bq 


(4.5) 


Because of its periodicity along z, the system can be paired to its Poincare map and 
therefore we will deal with a 2D discrete time system instead of a 3D continuous one. 
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Figure 5. Poincare plot of the magnetic field lines obtained taking the intersections between 
the Magnetic field lines and the snrface 2 = iz/2 with t — IISta- 


The resonant condition Beq ■ Vtf' = 0 leads to the definition of the resonant surfaces 
X = Xsi where the reconnection process takes place: 

d-ipeqix) d^ijdz 

dx dip^/dy kyi 

The initial perturbation considered by Borgogno et al. (2005) consists of two contributions 
with different wave number pairs (kzi, kyi)'. (1,0) and (1,1) respectively 

ip{x,y,z,t) = ipi{x,t) ejip{ikyiy + ikziz) + -ip 2 {x,t) exp{iky 2 y + ikz 2 z). (4.7) 

The amplitude of the perturbation 'ipi is of order 10“^ while ip 2 is of order 10“®. Magnetic 
islands are induced around resonant surfaces x = Xsi, i.e. Xsi = 0 and Xs2 = 0.71. When 
magnetic islands are sufficiently large the islands start interacting. Different modes of 
higher order are generated and the magnetic field topology becomes chaotic. The whole 
process can be visualized with the Poincare map technique which at any fixed t draws a 
snapshot of the magnetic field lines passing through a fixed z plane. Following (Borgogno 
et al. 2005) we consider the system after a time lapse of t = 415rA obtaining the Poincare 
map shown in Figure 5 with the intersections between the magnetic field lines and the 
z = Lz/2 surface. 

4.2. Numerical results and eomparison with FTLE ridges 
We look for the LCS of the system defined as the most repelling material lines. The 
finite transport barriers of the system also include the most repelling material fines 
obtained integrating backward in time, i.e. the attractive structures. We exploit the 
symmetry of the system under reflection over y = 0 to compute only the repulsive 
structures obtaining the attractive ones through their reflection. The algorithm used to 
integrate the dynamical system defined above is described by Onu et al. (2015). The 
result is shown in Figure 6 where red lines represent the repulsive structures while the 
blue ones the attractive ones. The domain can be split into different regions which have 
the LCS as boundaries. We expect fast mixing processes inside these regions. The number 
of material lines drawn depends on the filtering parameter chosen as explained by Onu 
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Figure 6. Plot of the most repelling material lines in red and of the most attractive ones in 
blue on the z = Z/z/2 plane superimposed to the Poincare map of the system. There is good 
agreement between the structures of the map and the material lines. 


et al. (2015). Lowering the filtering we obtain additional lines as depicted in Figure 7. A 
comparison between the structures obtained and the FTLE ridges shown in (Borgogno 
et al. 2011 a,&) shows a significant qualitative agreement. 


5. Conclusions 

This article is meant as an introduction to the study of Lagrangian Coherent Structures 
in chaotic magnetic field configurations. It aims to: 

1) stress the role that Lagrangian Coherent Structures can play in the description of 
transport phenomena in magnetically confined plasmas, 

2) introduce the plasma physics reader to a debate, occurring mostly in the fluid dynamics 
and oceanographic communities, on the proper operational definitions of these structures, 

3) exemplify the differences between the two definitions given by (Haller 2011) and by 
(Shadden et al. 2005) on a simple, analytically solvable, case, 

4) present a preliminary numerical comparison between the results given in previous 
works (Borgogno et al. 2011 a,&) where the first definition given byShadden et al. (2005) 
was used and a recalculation of the same structures based on the corrected definition 
given by Haller (2011). 

These results appear to be promising enough to start implementing a wider investi¬ 
gation of the applicability of LCS to magnetic configurations without the constraint of 
periodicity which is inherent to the Poincare map method. 
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Figure 7. As in Figure 6 but with a lower filtering parameter. Repelling structures are 
marked in red while the attracting ones with other colors. 


REFERENCES 

Avinash. K, Bulanov, SV, Esirkepov, T, Kaw, P, Pegoraro, F, Sasorov, PV & Sen, A 
1998 Forced magnetic field line reconnection in electron magnetohydrodynamics. Physics 
of Plasmas (1994-present) 5 (8), 2849-2860. 

Borgogno, D, Grasso, D, Pegoraro, F & Schep, TJ 2008 Stable and unstable 
invariant manifolds in a partially chaotic magnetic configuration generated by nonlinear 
reconnection. Physics of Plasmas (1994-present) 15 (10), 102308. 

Borgogno, D, Grasso, D, Pegoraro, F & Sghep, TJ 2011a Barriers in the transition to 
global chaos in collisionless magnetic reconnection, i. ridges of the finite time lyapunov 
exponent field. Physics of Plasmas (1994-present) 18 (10), 102307. 

Borgogno, D, Grasso, D, Pegoraro, F & Schep, TJ 20116 Barriers in the transition to 
global chaos in collisionless magnetic reconnection, ii. field line spectroscopy. Physics of 
Plasmas (1994-present) 18 (10), 102308. 

Borgogno, D, Grasso, D. Porcelli, F, Galifano, F, Pegoraro, F & Farina, D 2005 
Aspects of three-dimensional magnetic reconnection. Physics of Plasmas (1994-present) 
12 (3), 032309. 

Ghian, a G-L, Rempel, EL, Aulanier, G, Schmieder, B, Shadden, SG, Welsch, BT & 
Yeates, AR 2014 Detection of coherent structures in photospheric turbulent flows. The 
Astrophysical Journal 786 (1), 51. 

Goulliette, G, Lekien, F, Paduan, JD, Haller, G & Marsden, JE 2007 Optimal 
pollution mitigation in monterey bay based on coastal radar data and nonlinear dynamics. 
Environmental science & technology 41 (18), 6562-6572. 

Haller, G. 2011 A variational theory of hyperbolic lagrangian coherent structures. Physica D: 
Nonlinear Phenomena 240 (7), 574-598. 

Haller, G. 2015 Lagrangian coherent structures. Annual Review of Fluid Mechanics 47, 137- 
162. 





































14 


M. V. Falessi, F. Pegoraro and T.J. Schep 


Hamlington, P E, Poludnenko, AY & Oran, ES 2011 Interactions between turbulence and 
flames in premixed reacting flows. Physics of Fluids (1994-present) 23 (12), 125111. 

Huhn, F, Kameke, a, Perez-Munuzuri, V, Olascoaga, MJ & Beron-Vera, FJ 2012 The 
impact of advective transport by the south indian ocean countercurrent on the madagascar 
plankton bloom. Geophysical Research Letters 39 (6). 

Onu, K, Huhn, F & Haller, G 2015 Lcs tool: a computational platform for lagrangian coherent 
structures. Journal of Computational Science . 

Peacock, T & Haller, G 2013 Lagrangian coherent structures:: The hidden skeleton of fluid 
flows. Physics today 66 (2), 41. 

Peng, J & Dabiri, JO 2009 Transport of inertial particles by lagrangian coherent structures: 
application to predator-prey interaction in jellyfish feeding. Journal of Fluid Mechanics 
623, 75-84. 

PoRCELLi, F, Borgogno, D, Galifano, F, Grasso, D, Ottaviani, M & Pegoraro, F 2002 
Recent advances in collisionless magnetic reconnection. Plasma physics and controlled 
fusion 44 (12B), B389. 

Rappazzo, AF & Parker, EN 2013 Current sheets formation in tangled coronal magnetic 
fields. The Astrophysical Journal Letters 773 (1), L2. 

Rempel, el, Chian, A C-L, Brandenburg, A, Munoz, PR & Shadden, SC 2013 Coherent 
structures and the saturation of a nonlinear dynamo. Journal of Fluid Mechanics 729, 
309-329. 

Shadden, SC, Lekien, F & Marsden, Jerrold E 2005 Definition and properties of lagrangian 
coherent structures from finite-time lyapunov exponents in two-dimensional aperiodic 
flows. Physica D: Nonlinear Phenomena 212 (3), 271-304. 

Shadden, SC & Taylor, CA 2008 Characterization of coherent structures in the cardiovascular 
system. Annals of biomedical engineering 36 (7), 1152-1162. 

Tang, W, Mathur, M, Haller, G, Hahn, DC & Ruggiero, FH 2010 Lagrangian coherent 
structures near a subtropical jet stream. Journal of the Atmospheric Sciences 67 (7), 
2307-2319. 




